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Abstract. We present a comprehensive rheological study of a suspension of tfiermosensitive particles dispersed in water. 
Tiie volume fraction of tliese particles can be adjusted by the temperature of the system in a continuous fashion. Due to the 
finite polydispersity of the particles (standard deviation: 17%), crystallization is suppressed and no fluid-crystal transition 
intervenes. Hence, the moduli G' and G" in the linear viscoelastic regime as well as the flow curves (shear stress C7{f) as 
the function of the shear rate y) could be measured in the fluid region up to the vicinity of the glass transition. Moreover, 
flow curves could be obtained over a range of shear rates of 8 orders of magnitude while G' and G" could be measured 
spanning over 9 orders of magnitude. Special emphasis has been laid on precise measurements down to the smallest shear 
rates/frequencies. It is demonstrated that mode-coupling theory generalized in the integration through transients framework 
provides a full description of the flow curves as well as the viscoelastic behavior of concentrated suspensions with a single set 
of well-defined parameters. 
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INTRODUCTION 

The flow behavior of concentrated suspensions of hard sphere particles in a shear flow is a long-standing problem 
in colloid physics [Russel et al. (1989), Larson (1999)]. The most prominent feature of these systems is the strong 
increase of the zero-shear viscosity with increasing volume fraction ^. Concomitantly, there is a strong increase of the 
storage modulus G'. Moreover, with increasing shear rate a marked shear thinning sets in where the relative viscosity 
77/775 (775: viscosity of solvent) is lowered considerably. These general features are also observed for solutions of soft 
colloidal particles as e.g. a solution of star polymers in a suitable solvent [Helgeson et al. (2007)]. Suspensions of hard 
spheres or of particles characterized by a sufficiently steep mutual repulsion may hence be treated as the most simple 
model system for the rheological study of the flow behavior of colloidal systems. 

By now, it is well understood that the dynamics of quiescent suspensions of hard spheres can be understood 
and modeled quantitatively in terms of the mode-coupling theory [van Megen and Pusey (1991), van Megen and 
Underwood (1994)]. Mode-coupling theory [Gotze and Sjogren (1992)] was shown to explain quantitatively the 
structural arrest that takes place in colloidal dispersions above the volume fraction (j)g w 0.58. Moreover, this value 
of (j)g has turned out to be in full agreement with experimental studies of the zero-shear viscosity 770 [Meeker et al. 
( 1 997), Cheng etal. (2002)] . Here it has been demonstrated that the strong upturn of 770 occurs much below the volume 
fraction of 0.63 that would refer to random close packing of the spheres. First theoretical calculations of 770 within the 
frame of MCT fully supported this view [Banchio et al. (1999)]. 

Fuchs and Cates (2002, 2003) recently demonstrated that MCT can be generalized to treat the dynamics of colloidal 
suspension in shear flow. They developed the integrations throught transients (ITT) approach, which describes the 
competition between enforced deformation or flow and intrinsic structural relaxation. Thus, MCT can be used to 
model the distortion of the structure of suspensions by a flow field at finite shear rate 7. One central point of this theory 
is the melting of the glass by shear: flow leads to a glass-to-fluid transition in the suspension and a unique stationary 
state can be reached again despite the fact that the volume fractions of the particles in the suspension are above It 
is presumed that this holds as long as local 'free volume' is available below random close packing. Flow curves giving 



FIGURE 1. Schematic representation of the thermosensitive core-shell particles used in this study: The polymer chains are affixed 
to the surface of the core. Immersed in cold water the thermosensitive network will swell. The degree of swelling can be reduced 
by raising the temperature. At 32 "C a volume transition takes place in which most of the water is expelled from the network. 



the relative viscosity rj/ris as function of shear rate are predicted to present well-defined states that are independent 
of the previous history of the sample. In addition to these predictions from ITT-MCT pertaining to the non-linear 
flow of suspensions, MCT is capable of describing the linear-viscoelastic behavior of concentrated dispersions. In 
this way MCT offers a complete modeling of the rheological properties of dispersions at high packing fractions (j) or 
for strong particle interactions close to vitrification. First comparisons with simulations indeed demonstrated that the 
theory captures all of the salient features of the flow behavior of flowing suspensions [Varnik and Henrich (2006)]. 

Recently, it has been demonstrated that suspensions of thermosensitive particles present excellent model systems 
for studying the flow of concentration suspensions for a quantitative comparison of the theory of Fuchs and Gates 
to experimental data [Fuchs and Ballauff (2005), Crassous et al. (2006a, 2008a)]. The particles consist of a solid 
core of polystyrene onto which a thermosensitive network of poly(N-isopropylacrylamide) (PNIPAM) is attached 
([Crassous et al. (2006b, 2008b)], see Figure[T]i. The PNIPAM shell of these particles swells when immersed in cold 
water (below 32 "C). Water gets expelled at higher temperatures leading to a considerable shrinking. Thus, for a given 
number density the effective volume fraction ^^ff can be adjusted within wide limits by adjusting the temperature. 
Senff et al. (1999) were the first to demonstrate the use of these particles as model system for studying the dynamics 
in concentrated suspensions [Senff et al. (1999)]. Similar microgel particles have been used to study the dynamics 
of concentrated suspensions and the aging in these systems [Purnomo et al. (2006, 2007)]. The advantage of these 
systems over the classical hard sphere systems are obvious: Dense suspensions can be generated in situ without shear 
and mechanical deformation. Moreover, all previous history of the sample can be erased by raising the temperature 
and thus lowering the volume fraction to the fluid regime. Hence, subtle effects as e.g. aging can be studied by these 
systems in a nearly ideal fashion [Purnomo et al. (2006, 2007)]. Moreover, microgels can be used as models for pastes 
and gels in which the effective volume fraction of the particles exceeds unity by far [Seth et al. (2006)]. 

Previous investigations using dispersions of thermosensitive core-shell particles have demonstrated that mode- 
coupling theory is indeed able to model the flow curves as well as the linear viscoelasticity of concentrated suspensions 
[Crassous et al. (2006a, 2008a)]. In particular, this theory allows us to calculate G' and G" as function of the frequency 
(£) from data that have been obtained from the analysis of flow curves [Crassous et al. (2008a)]. However, a problem 
that restricted this comparison was crystallization due to the narrow size distribution of the particles. Crystallization 
presents a major obstacle since the rate of nucleation is rather fast for volume fractions even far below the volume 
fraction of glass transition (0^ « 0.58). The crystalline phase could be removed by shear melting and the fluid state 
could be super-cooled for a given time. The time window of nucleation, however, limited the experimental studies to 
frequencies and/or shear rates not too small [Crassous et al. (2006a, 2008a)]. A fully conclusive comparison of MCT 
with experimental data, however, requires that measurements extend down to the smallest frequencies and shear rates 
[Fuchs and Baflauff (2005), Crassous et al. (2008a)]. 

Here we present a full comparison of theory and experiment using a thermosensitive suspension having a slightly 
broader size distribution that prevents crystallization. Hence, all measurements can now be extended to very low shear 
rates or frequencies that may take hours without any disturbance by crystallization. It will be demonstrated that MCT 
provides indeed a full description of the rheological behavior of suspension under homogeneous flow. Special attention 
wiU be paid to shear banding [Caflaghan (2008), Dhont and Briels (2008), ManneviUe (2008), Olmsted (2008)] and 



other effects that may come into play for concentrated suspensions [Isa et al. (2007)J and disturb the quantitative 
comparison of theory and experiment. 

The paper is organized as follows: In the following section we give a brief survey of the main predictions of MCT. 
In particular, the equations used for the evaluation of data as well as the fit parameters necessary for a comparison 
of MCT with experimental data are discussed. For a extensive description of MCT the reader is deferred to a recent 
paper [Crassous et al. (2008a)]. After a section devoted to the description of the systems and the measurements, a 
comprehensive comparison of MCT with flow curves as well as with linear viscoelastic properties will given. A major 
focus is the low frequency and long time behavior which had not been accessible before. Special attention will be paid 
to the fact that MCT is capable of describing both the flow curves as well as G' and G" as functions of frequency. A 
brief Conclusion wiU sum up the main findings. 



MODE-COUPLING THEORY OF FLOWING SUSPENSIONS: MAIN PREDICTIONS 



The central points and main predictions of MCT and its generalization to flowing dispersions in the integration through 
transients (ITT) approach have recently been discussed in detail [Crassous et al. (2008a)]. Hence, it will suffice to 
discuss the points that will be used in the analysis of the experiments. Structural relaxation, which is the central topic 
of MCT, leads to internal relaxation processes with an intrinsic relaxation time T that far exceeds the dilute limit 
Brownian diffusion time. The ITT-MCT approach to sheared fluids considers the effect of shear when the dressed 
Peclet or Weissenberg number Pe = jx becomes appreciable, while the bare Peclet number Pe^ = 6nrisRjj'Y/ {ksT) 
formed with the solvent viscosity t/^ and the hydrodynamic particle radius Rh should be much lower than unity for the 
approach to hold. 

The main point of the MCT of quiescent dispersions is the slowing down of the dynamics at the glass transition. 
Criticid values for the particle concentration, the temperature or other thermodynamic control variables exist, where 
a glass transition separates fluid from yielding glassy states. In fluid states, a first Newtonian region is reached if the 
shear rate or the frequency is low enough; the Newtonian viscosity JTo = '7(7^ 0) is finite. Non-Newtonian flow, 
that is shear thinning, is observed at higher shear rates corresponding to a finite Weissenberg number Pe « 1 . MCT 
recovers Maxwell's relation Tjo <x t, so that the criterion Pe « 1 translates into T7o7» GL, where = G'{(0—^ 0) is the 
(zero-frequency) shear modulus at the glass transition. For glassy states at higher densities or interaction parameters, 
Newtonian flow is not observed anymore but the shear stress a is predicted to approach a plateau value for sufficiently 
low shear rates y; a dynamic yield stress follows (7+ = (7(7—+ 0) > 0. It corresponds to an infinite Newtonian viscosity. 

A central point of MCT as applied to flowing suspensions [Fuchs and Cates (2002)] is shear melting: Finite shear 
rates will lead to a melting of the glass into a fluid. As a consequence of this, flow curves given by the shear stress a 
as function of the shear rate 7 present well-defined states for a given volume fraction even above the glass transition. 
Hence, a (7) is independent of the previous history of the sample, under the condition that the stationary long time 
limit has been achieved. Previous experimental studies strongly suggest that this central prediction of theory is found 
indeed [Crassous et al. (2006a, 2008a)]. However, in all previous experiments crystallization intervened at sufficiently 
low shear rates and this point is in need of further experimental elucidation. 

In the following, the main equations used for the evaluation of the rheology data within MCT are briefly summed 
(y) 

up. Within the schematic -model [Fuchs and Cates (2003)] Brownian motion of the particles is captured in terms 
of a density correlator which is normalized/initiated according to ^ 0) = 1 — R. MCT relates this correlation 
to the memory function m{t) by a memory equation 



The parameter F describes the microscopic dynamics of the density correlator and will depend on structural and hy- 
drodynamic correlations. The memory function describes stress fluctuations which become slower when approaching 

(y) 

the volume fraction of the glass transition. In the F^l -model the memory function is written as 




(1) 



m{t) = 



Vi<i>{t) + V2'S>^ 



it) 



(2) 



1 + (7f/7c)' 



This model, for the quiescent case 7 = 0, had been suggested by Gotze in 1984, [Gotze (1984, 1991)] and describes the 
development of slow structural relaxation upon approaching a glass transition; it is reached by increasing the coupling 



vertices/ parameters vi and V2 to their critical values Vj ~ (2A — 1)/^^ and = where the so-called exponent 

parameter A with 1/2 < A < 1 parametrizes the (type-B) glass transition line. Under shear an explicit time dependence 
of the couplings in m(t) captures the accelerated loss of memory by shear advection [Fuchs and Gates (2002, 2003, 
2008)] . Shearing causes the dynamics to decay for long times, because fluctuations are advected to smaller wavelengths 
where small scale Brownian motion relaxes them. The parameter 7^ sets the scale for the accumulated strain yt to 
matter 

Equations ([T] |2]i will be used with the choice of vertices V2 = = 2 (corresponding to A = l/\/2), and vi = 
Vj + £ (1 — f^)/ , where Vj = 0.6258. The critical glass form factor f^, which takes the value = (1 — A) = 0.293 
for this choice, describes the frozen-in structure obtained from <I>(f °°,e =0,7=0) = /' at the glass transition in 
the absence of shear. The glass transition lies at e = 0, where the Newtonian viscosity diverges. Under shear, glass 
states are molten, but the transition still separates shear thinning fluid states, e < 0, where a Newtonian viscosity exists, 
from yielding glassy states e > 0, where a dynamic yield stress exists. The specific choice of this transition point is 
motivated by (i) its repeated use in the literature, and (//) because it predicts a divergence of the viscosity at the glass 
transition, whose power is close to the one of the microscopic MCT calculation for short ranged repulsions: 

770 oc T oc (-e)-l' , with 7 = 7(A = 1 /^/2) = 2.34 . 

Within microscopic MCT for hard spheres using the Percus-Yevick structure factor 7 = 2.46 was found [Gotze and 
Sjogren (1992)]. 

The correlator <I>(f) obtained by numerical solution of Eq. [T]may then be used to calculate the generalized shear 
modulus g{t,Y) according to 

g{t,j) = Va^\t) + ^^5{t~Q+). (3) 

The parameter Tjoo characterizes a short-time high frequency viscosity and models viscous processes which require no 
structural relaxation. Together with F, it is the only model parameter affected by solvent mediated interactions. Steady 
state shear stress under constant shearing, and viscosity then follow via integrating up the generalized modulus: 

CT = 77 7= 7 / dtg{t,y) = y [ dt Va^^it) + yri„ . (4) 
Jo Jo 

When setting the shear rate 7 = in Eq.|2]so that the schematic correlator belongs to the quiescent equilibrium system, 
the frequency dependent moduli of the linear response framework follow from Fourier transforming: 

G'{co) + iG"{co) = ico [ dt e'''^' g{t,Y = 0) =ico [ dt e'''^' Va ^^it)\^n + icori„ (5) 
Jo Jo 

Because of the vanishing of the Fourier-integral in Eq. |5] for high frequencies, the parameter rjoc can be directly 
identified as high frequency viscosity: 

Um G"(a))/« = 77f , with 77^ = t]^ . (6) 

(B— ►oo 

At high shear, on the other hand, Eq. [2] leads to a vanishing of m(f), and Eq.[T] gives an exponential decay of the 
transient correlator, <I>(f) e^^' for 7^0. The high shear viscosity thus becomes 

77i=lim(T(7)/7=77o„ + ^ = C + ^. (7) 

Eq. I2] connects the high-shear viscosity with the high-frequency limit of the viscosity in an analytical manner 

Moreover, both rjl and rj^ are accessible from experimental data and Eq.|7]provides an additional check for internal 

(y) 

consistency of the fit within the frame of the schematic model. Representative solutions of the 2 -model have been 
discussed recently [Crassous et al. (2008a), Fuchs and Gates (2003), Hajnal and Fuchs (2008)]. 



EXPERIMENTAL 



The particles consist of a solid core of poly(styrene) onto which a network of crosslinked poly(N-isopropylacrylamide) 
(PNIPAM) is affixed [Grassous et al. (2006b, 2008b)]. The degree of crosslinking of the PNIPAM shell effected by 



the crosslinker N,N'-methylenebisacrylamide (BIS) was 2.5 Mol %. The core-shell type PS-PNIPAM particles were 
synthesized, purified and characterized as described in [Dingenouts et al. (1998)]. Immersed in water the shell of 
these particles is swollen at low temperatures. Raising the temperature above 32 "C leads to a volume transition within 
the shell. Cryogenic transmission electron microscopy and dynamic light scattering have been used to investigate the 
structure and swelling of the particles [Crassous et al. (2006b, 2008b)]. Further experimental details on cryogenic 
transmission electron microscopy are given in [Crassous et al. (2006b)]. 

The average size and size distribution of the particles was determined with the CPS Disc Centrifuge from LOT, model 
DC 24000 in a aqueous glucose density gradient between 3 and lwt% as medium at room temperature. The number- 
average diameter determined at 25 "C is 131 nm and the standard deviation of the size distribution is 17 %. 
Dynamic light scattering (DLS) was done with a Peters ALV 800 light scattering goniometer The temperature was 
controlled with an accuracy of 0.1 "C. In the range of temperatures used for the measurements, Rh could be linearly 
extrapolated between 6 and 25"C {Rh = -0.7796 T H-102.4096 with T the temperature in "C) (see [Crassous et al. 
(2006a)]). Fig.|2j5 displays Rh as the function of temperature. Moreover, this plot demonstrates that Rh can obtained 
very precisely. Since the subsequent reduction of the data (see below) involves a scaling of the data by Rj^, this fact is 
of great importance. 

Density measurements have been made using a DMA-60 densiometer supplied by Paar (Graz, Austria). The 
density for the core-shell particles were determined as follows: p = l.llSA-gcm^^ at 25"C, 1 .1239 g cm^^ at 20"C, 
1.124-6gcm^^, at 18"C and 1.1259gcm^^ at IT'C. These densities can be described with the function p = 6.800- 
10^^ • + 1.129 with the temperature T in "C. In all experiments reported herein S.lO^^molL^^ KCl were added 
to the suspensions in order to screen the remaining electrostatic interactions caused by a few charges attached to the 
surface of the core particles. Previous studies have shown that this screening is fully sufficient. Thus, within the range 
of volume fractions under consideration here the thermosensitive particles can be regarded as hard spheres (see the 
discussion of this point in [Crassous et al. (2006a)]). 

The flow behavior and the linear viscoelastic properties for the range of the low frequencies were measured with 
a stress-controlled rotational rheometer MCR 301 (Anton Paar), equipped with a Couette system (cup diameter: 
28.929mm, bob diameter: 26.673mm, bob length: 39.997mm). For comparison the double gap geometry (inner cup 
diameter: 23.820mm, inner bob diameter: 24.660mm, outer bob diameter: 26.668 mm, outer cup diameter: 27.590mm 
and bob length: 40.00 mm) was used as well. Measurements have been performed on 12 mL (Couette single gap 
geometry) and 4mL (Couette double gap geometry) solution and the temperature was set with an accuracy of ± 
0.05 "C. Deformation tests were done from a deformation of 0.1 to 100% in the case of the dynamic deformation 
test with a logharitmic time ramp of 200 to 20 s at IHz. The shear stress a versus the shear rate 7 (flow curve) was 
measured after a pre-shearing of 7= 100 i^' for two minutes and a waiting time of lOi. The flow curves were obtained 
setting 7 to a given value, first with increasing 7 from 7 = 10^^ — 10^ s^^ with a logarithmic time ramp from 1500 to 
20 s. Then after a waiting time of 10.? 7 was decreased again. 

The frequency dependence of the loss G" and elastic G' moduli has been measured after a two minute pre-shear 
at a shear rate of 100.?^' and 10.? waiting time for 1% strain from 10 to 10^'^ Hz with a logarithmic time ramp 
from 20 to 1500i. The dependence upon the strain has been checked and confirmed that all the measurements were 
performed in the Unear regime. G' and G" at elevated frequencies were obtained using the Piezoelectric vibrator 
(PAV) [Crassous et al. (2005)] and cylindrical torsional resonator [Deike et al. (2001), Fritz et al. (2003)] supplied by 
the Institut fiir dynamische Materialpriifung, Ulm, Germany. The PAV was operated from 10 to 3000 Hz. The solution 
is placed between two thick stainless steel plates. One remains static whereas the second is cemented to piezoelectric 
elements. The gap was adjusted with a 50 m ring. A set of piezoelectric elements is driven by an ac voltage to induce 
the squeezing of the material between the two plates. The second set gives the output voltage. Experimental details 
concerning this instrument are given elsewhere [Crassous et al. (2005)]. The cylindrical torsional resonator (TR) 
used was operated at a single frequency (26 kHz). The experimental procedure and the evaluation of data have been 
described recently [Deike et al. (2001), Fritz et al. (2003)]. 
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FIGURE 2. a) Cryogenic Transmission Electron Microscopy micrograph of the core-shell particles in water, b) Dynamic light 
scattering measurements of the latex at varius temperatures. Filled squares represent the hydrodynamic radius measured in aqueous 
QmmolL-^ KCl, hollow squares measured in water. All rheological experiments presented here have been done using particles 
dissolved in aqueous 0.05 mo/ L^' KCl solution [Crassous et al. (2006a, 2008a)]. In this way possible electrostatic interaction 
between the particles due to the residual charges affixed to the core particles can be screened totally. Below the volume transition 
the decrease of the hydrodynamic radius Rh giving the outer boundary of the network is linear in very good approximation {Rh = 
-0.7796 -T +102.4096 with T being the temperature in "C; see solid line in graph). 



RESULTS AND DISCUSSION 
Characterization 

Recent work has shown that Cryo-TEM is a excellent tool for the analysis of the core-shell particles used for the 
present study [Crassous et al. (2006b, 2008b)] because this method allows us to visualize the particles in-situ, that is, 
directly in the aqueous suspension. Fig.|2]displays a micrograph of the suspension obtained by cryo-TEM. The core- 
shell structure of the particles is clearly visible as akeady demonstrated by previous work [Crassous et al. (2006b, 
2008b)]. Since the polydispersity of the present suspension is considerably higher than the one of the dispersion used 
in our previous study [Crassous et al. (2008a)], the system does not crystallize anymore even if kept for a prolonged 
time at volume fractions above 50%. However, Fig.|2]demonstrates that the suspension used here is still a well-defined 
model system. Hence, it can be used for a comprehensive comparison of experimental data with MCT. Attempts to 
suppress crystallization by mixing core-shell particles of slightly different size failed to produce reliable results since 
crystallization in these systems was found to be only slowed down but not suppressed totally. 

The effective volume fraction i^^ff could not be calculated as for monodisperse particles by precise measurements 
of the hydrodynamic radius Rh, the core radius and the mass ratio of core and shell as done in [Crassous et al. 
(2008a)]. Unfortunately, polydispersity leads to an error which is too large and no meaningful data could be obtained. 
Furthermore, viscosity measurements in the dilute state and the use of the Batchelor-Einstein equation for the relative 
zero shear viscosity could not be applied because of the appreciable polydispersity, too. Therefore the effective volume 
fraction was obtained by comparing the relative viscosities at infinite shear rate with the monodisperse core-shell 
system used in [Crassous et al. (2008a)]. This will be discussed in further detail in conjunction with Fig.|6]below. 



Reversibility of flow curves 

A central prediction of MCT is the full reversibility of the flow curves, that is, of the shear stress a as the function of 
the shear rate y. This is evident for any volume fraction below the volume fraction <^g of the glass transition. Here the 
system is ergodic and a well-defined steady state must be reached for sufficiently long time. This point is well borne 
out of all experiments done so far [Crassous et al. (2006a, 2008a)]. However, a totally different situation arises above 
^g. In the glassy state, the suspension becomes non-ergodic and the flow curves may become ill-defined states. Hence, 
for concentrated suspensions of large non-Brownian particles dispersed in a highly viscous medium, the measured 
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FIGURE 3. a) Reversibility of the flowcurves for a 7.87 w/% suspension at 12°C with two different geometries: Couette single 
gap - geometry (crosses) and double gap - geometry (squares), b) Reversibilty of the flowcurves by comparing jump experiments 
and the flowcurve (solid line) measured from higher to lower shear rates at SA2wt% and 14.75°C. The steady state shear stress 
values of jumps from different shear rates to one as shown by the arrows in the graph are plotted for two examples: jump from 
lOOi-"' (full circle) to 0.000127 5^' (full square) and jump from O.OIO85"' (hollow circle) to 0.000127 5"' (full square), b) Inset: 
the same jumps are plotted time-dependent. The different jumps are plotted like followed: jump from lOOi^' to 0.000127^^' (thick 
solid line) and jump from 0.0108^^' to 0.000127 (thick dashed line). The thin short dashed line in the inset represents the shear 
stress value at 0.000127 .s^' of the flowcurve for the same concentration and temperature. 



shear stress is found to be dependent on the previous history of the sample [Heymann and Aksel (2007)]. MCT, 
however, predicts that shear melting sets in at any finite shear rate 7. Therefore, shear flow restores the ergodicity in 
the system and thermal motion of the spheres can erase all traces of the previous deformations of the sample if flow 
is applied over a sufficiently long time interval. Two types of experiments need to be done to explore these basic 
principles in further detail: 

i) Each data point giving (7(7) must be independent regardless of the shear history. Thus, flow curves (7(7) should be 
fully reversible regardless whether the data are taken by increasing rate or decreasing the shear rate. The reversibility 
of the flow curves has already been demonstrated in previous work. However, the suspensions analyzed previously 
[Crassous et al. (2006a, 2008a)] began to crystallize during the long measurements needed for exploring the dynamics 
at smallest shear rates 7. To overcome this problem, the present system is sufficiently polydisperse, crystallization 
is suppressed. Fig. |3] displays the resulting flow curves down to Pcq = 10^^. The onset of a plateau is clearly seen. 
Additional experiments with long waiting times demonstrated that crystallization does not occur in the present system. 
Moreover, Fig. points to another important fact: (7(7) is entirely independent of the geometry of the rheometer 
used for studying the flow behavior Fig. [3^ demonstrates this fact by overlaying the resulting ff (7) obtained from 
a single gap and double gap Couette geometry. If shear banding would play any role [Callaghan (2008), Dhont and 
Briels (2008), Manneville (2008), Olmsted (2008)], the resulting a{y) should be different and depend definitely on the 
width of the gap of the instrument. The difference in curvature for the two experiment geometries would have resulted 
in different onsets of shear banding. Also, shear banding in the double-gap Couette system would have resulted in a 
complex interplay of instabilities in the two gaps. None of this is observed here. 

ii) The full reversibility of the flow curves are demonstrated also by a second type of experiment (Fig.|3}'): Here a 
stationary state is first obtained for a given shear rate 7. Then a much lower shear rate is adjusted and the resulting 
shear stress is monitored as a function of time. Fig. shows that (7(7) converges to the stationary value measured 
previously (solid line in Fig.fSj)). The inset displays the evolution of (7(7) with time. Evidently, a long time may be 
needed until the stationary value is reached. However, all data demonstrate that the full reversibility of (7(7) is fully 
supported by our experiment also above the glass transition. One of the central predictions is thus fully corroborated. 
Deviations observed previously [Crassous et al. (2008a)] were solely due to crystallization of the sample which is 
suppressed in the present system. In this way flow curves may to a certain extent be compared to thermodynamic 




FIGURE 4. Comparison of the dynamic strain sweep at IHz (C: filled squares, G" open squares, a solid line) and a continuous 
strain sweep with two different measurement times for each point (20s per point: triangles down and 2005 per point: circles) for a 
8.90w/% suspension at 17°C (0^// = 0.627). 



results, they are fully independent of the path leading to a given state. 

A point to be considered further is the transition from the linear to the non-linear viscoelastic regime in the glassy 
regime. Fig. |4] displays G', G" as well as a for different deformations at a rate of 1 Hz. The data show that the linear 
regime is preserved up to a strain of nearly 10%. Moreover, G'y measured for different strains 7 coincides directly 
with the measured stress a. Hence, within this range of 7, the glassy suspensions behaves as a homogeneous solid as 
expected. Again, shear banding or possible inhomogeneities within the gap should disturb this type of measurement. 
Fig.|4]demonstrates that this is not seen. Obviously, these experiments are not fully sufficient to rule out shear banding 
entirely. However, the good agreement between two entirely different types of experiments with theory to be discussed 
farther below points clearly to the fact that shear banding does not disturb the measurements reported here. 

Beyond a strain of ca. 10 %, shear melting of the glass starts and the measured stress cj is considerably lower than 
the value calculated from G'. Fig. |4] indicates that the linear viscoelastic regime is extended over a sufficiently large 
range of deformations. 



Comparison with mode coupling theory: Schematic model 

A thorough comparison of ITT-MCT to experimental data must be related not only to the flow curves discussed in the 
previous section but also to data obtained for the linear viscoelasticity, that is for G' and G" as function of the frequency 
CO. A previous comparison [Crassous et al. (2008a)] has demonstrated that both sets of data can indeed be described by 
a generalized modulus ^(f , 7) as introduced in the section Theory. However, at smallest Peo crystallization intervened 
and the comparison was restricted to intermediate to high shear rates/frequencies. Crystallization is suppressed for 
the present system and a fully quantitative comparison of theory and experiment now becomes possible down to 
Peo « 10^^ in case of the flow curves. G' and G" can be obtained down to Pea « 10^^. Here, the frequency dependent 
Peclet (or Deborrah) number, defined as Pe^) = 67rrisRjf(o/{kBT), compares the frequency with the time it takes a 
particle to diffuse its radius Rh at infinite dilution. At high densities, structural relaxation shows up at Pcoj ^ 1 . 

In the following we shall compare the experimental flow curves and G' and G" with theory for a given effective 
volume fraction ^eff- Fig.|5]gives this comparison for five different effective volume fractions. On the left-hand side 
the flow curves (1(7) are presented as function of the bare Peclet number Peo, on the right-hand side G' and G" are 
displayed as the function of Pe^, now calculated with the respective frequencies co. Table [T] gathers the respective 
volume fractions together with the other experimental conditions. Note that G' and G" have been obtained over nearly 
7 orders of magnitude in frequency by combining a Couette rheometer (single gap, double gap) with the PAV and the 
torsional resonator. 

The generalized shear modulus g{t,Y) (cf. Eq. |3]l presents the central theoretical quantity used in this fit. Within 
the schematic model, the vertex prefactor is kept as a shear-independent quantity. It can easily be obtained 
from the stress and modulus magnitudes. Hydrodynamic interactions enter through Tjoo, which can be obtained from 
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FIGURE 5. Left column: Reduced flow curves for different volume fractions (graphs in left column a, c, e, g and /). Hollow 
squares represent the '"up"' curve from low to high shear rates, filled squares represent the "'down"'curve from high to low shear 
rates. The solid lines are the results of the schematic model, the dashed line represent the pseudo power law behavior (see Table 2). 
Right column: Reduced frequency dependent moduli for different volume fractions (graphs b, d, f, h and j 
in right column). Full symbols/solid lines represent G', hollow symbols/dashed lines represent G". Squares 
mark the data points obtained by the Physica MCR 301, circles with the PAV, triangles up with the TR, 
the thick lines are the results of the schematic model, the thin lines the results of the microscopic model. 
Graphs in one row represent the continuous and dynamic measurements at one volume fraction, a and b were measured at 
8.90wr% and 25 °C (peff = 0.530 , c and d were measured at 7.87 w/% and 18°C (circles) ^^ff = 0.595 , e and/ at 8.90w/% and 
20°C <l)eff = 0.616 , g and A at 8.90wf% and 19°C 0^// = 0.625 , and i and ; at 8.90wf% and 17°C 0^// = 0.627 . 




measurements done at high frequencies, and through F, which can be obtained via Eq.|7]from measurements done at 
high shear rates. Given these three parameters, both the shapes of the flow curves as weU as the shapes of the moduH 
G' and G" may be obtained as function of the two parameters e and y/Yc- The former sets the separation to the glass 
transition and thus (especially) the longest relaxation time, while the latter tunes the effect of the shear flow on the 
flow curve. The characteristic changes of the shapes of the steady stress and the moduli are dominated by £ and y/Yc- 

The fitting of the data was done as follows: All measured quantities, namely ff, G', and G" were converted to 
the respective dimensionless quantities by multiplication with Rj^/ksT where Rh is the hydrodynamic radius at the 
respective temperature. Previous work has demonstrated that the effective volume fraction deriving from Rh is the 
only decisive thermodynamic parameter [Crassous et al. (2006a, 2008a)]. As already discussed above, 7 and (O are 
converted by ^nrfsR^i /ksT to the respective Peclet numbers. 

Evidently, both the reduced moduli, the Pe number, and the packing fraction depend on the effective particle 
volume R^fj. Therefore, this quantity must be measured with utmost precision. In our poly disperse sample, actually 
a distribution of values exists, whose variance we determined by disc centrifugation at low concentration. It is 
important to note, however, that the size distribution does not change with temperature. Moreover, in our experiments 
close to the glass transition, we may assume the size distribution to be (almost) density independent. Therefore, we 
have R]j{T) as single experimental control parameter, whose small change upon varying temperature T drives the 
system through the glass transition (see also the discussion of Fig. 6). 

Subsequently, the fit procedure is started at the smallest volume fraction and a trial value for the separation parameter 
e is chosen. This fit can be done very accurately since the resulting reduced moduli depend very sensitively on e. Once 
the overall shape of both the reduced stress a/?^ /ksT and of the reduced storage and loss modulus has been reproduced 
by a proper choice of e, the vertex parameter is obtained. For lower volume fraction the crossing point of G' and 



TABLE 1. Experimental parameters (concentration c, temperature T and calculated ^eff) of the measurements shown 

in Fig. [5] The fit parameter va, F, jc and r\oo of the Fj^ model as well as the calculated values r\L, (Tq and r\^) are written in 
terms of the respective normalization. 
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* calculated by Eq.|2] 

^ extrapolated from the fits using the ^['2 ' -model 

** directly measured for the three lower volume fractions, extrapolated for the two highest volume fraction from the fits using the 
-model 



G" may be used, near the glassy state the plateau value of o and G' can be used. Note that in all steps both sets of data, 
that is, the flow curves as well as G' and G" must be fitted with the same Vc- Next Y is determined from the dynamic 
data which are fully described with this choice of parameters. Given 7, the parameter jc (see Eq.|2]) is chosen so that 
the experimental flow curves are matched by theory. Then Tjo, ctq and r\L can be calculated. The results of these fits 
are gathered in Table 1 . We reiterate that the final set of parameters must lead to a description of both the flow curves 
as well as of the dynamic data. This requirement stresses the connection of the equilibrium structural relaxation at 
the glass transition to the non-linear rheology under strong flow, which is one of the central points of the ITT-MCT 
approach. 

Fig. |5] demonstrates that the rheological behavior of a non-crystallizing colloidal can be modeled in a highly 
satisfactory manner by five parameters that display only a weak dependence on the effective volume fraction of the 
particles. Increasing the effective packing fraction drives the system toward the glass transition, viz. e increases with 
^eff- Stress magnitudes (measured by Va) also increase with as do high frequency and high shear viscosities; 
their difference determines Y. The strain scale 7 remains around the reasonable value 10%. In spite of the smooth 
and small changes of the model parameters, the -model achieves to capture the qualitative change of the linear and 
non-linear rheology. The measured Newtonian viscosity increases by a factor around 10^. The elastic modulus G' at 
low frequencies is utterly negligible at low densities, while it takes a rather constant value around lOksT /E?[j at high 
densities. An analogous observation holds for the steady state shear stress o{j), which at high densities takes values 
around 0.3kBT /Rjj when measured at lowest shear rates. For lower densities, shear rates larger by a factor around 10^ 
would be required to obtain such high stress values. 

At volume fractions around 0.5 the suspension is Newtonian at small Peo (see Fig.|5]l. Evidently, these flow curves 
are fully reversible since the suspension is still in the fluid regime. The schematic model provides a perfect fit. 
Moreover, theory demonstrates in full agreement with experiment that the high-shear viscosities measured in shear 
flow and the high-frequency viscosities measured in the linear viscoelastic regime do not agree. As is evident from Eq. 
I2I both rjL and rj^ are related by Eq.|7]and the Cox-Merz rule does not hold for suspensions as expected. Approaching 
the glass transition leads to a characteristic S-shape of the flow curves and the Newtonian region becomes more and 
more restricted to the region of smallest Peo. Concomitantly, a pronounced minimum in G" starts to develop, separating 
the slow structural relaxation process from faster, rather density independent motions, while G' exhibits a more and 
more pronounced plateau. At the highest densities (Fig. |5]panels i and j), the theory would conclude that a yielding 
glass is formed, which exhibits a finite elastic shear modulus (elastic constant) Geo = G'{(0 ^ 0), and a finite dynamic 
yield stress, (7+ = (7(7-^ 0). The experiment shows, however, that small deviations from this glass like response exist 
at very small frequencies and strain rates. Description of this ultra-slow process requires extensions of the present 
ITT-MCT [Crassous et al. (2008a)], where additional physical mechanisms for stress relaxation need to be considered 
that are neglected here. 



100 

50 

V 20 
^ 10 

V 

5 
2 

0.40 0.45 0.50 0.55 0.60 0.65 

iff 

FIGURE 6. Relative viscosities in dependence of the effective volume fraction for a monodisperse latex [Crassous et al. (2008a)] 
(circles) and the polydisperse latex (this work; squares). The hollow symbols show the high frequency viscosities, the filled symbols 
the high shear viscosities. The relative high frequency viscosity r]^ can be described by Lionberger-Russel - equation (see Eq. 
4.22 in Lionberger and Russel (1994)) with a reduced effective volume fraction of 11.95%. Assuming the same relation for the 
polydisperse latex under consideration here, Tj^ was used to determine the effective volume fractions. The (l)gff obtained by this 
way are plotted in the inset versus the concentration c [g • niL^^] times the cubed hydrodynamic radius Rj^ [m^]. The solid line 
represents the linear scaling of effective volume fraction and c ■ Rj^ with a slope of 9.77- 10^ g and a intercept of zero. As a further 
check, the high-shear viscosities are plotted as well as a function of (peff- The filled circles show this quantity for the narrowly 
distributed latex [Crassous et al. (2008a)] whereas the filled squares denote the repsective quantities for the present system. Good 
agreement is seen thus corroborating the present way of determining <l)eff- 

Power laws 

Often power-law shear thinning is discussed, when the viscosity is plotted as function of shear rate instead of 
the flow curves, viz. stress versus shear rate. For fluid states, the double logarithmic presentation in Fig. |5] shows 
characteristically S-shaped flow curves, which have a point of inflection. Clearly, this flow curve being a more sensitive 
presentation as the usual Tj(7)-plots rules out the existence of a true power law behavior, which should show up as 
a straight line in the double logarithmic plot. Yet, at the point of inflection of the S-curve, and closely around it, a 
pseudo-power law may exist. Here rj (7) ^ jiP^^)^ where p is the slope of the tangent in the flow curve at the inflection 
point, p — dlna /dlnj- These tangents are indicated by dashed straight lines in the panels of Fig. 5. The pseudo- 
power law that corresponds to the one in the viscosity can describe the flow curve over more than two decades in shear 
rate. The exponents p (see Table |2]i change with density, and can be computed within the schematic model [Hajnal 
and Fuchs (2008)]. Nevertheless, Fig.|5]shows clearly that the description by power laws is by far inferior to the one 
provided by MCT. 



Effective volume fraction 

As akeady lined out in the section Experimental, the higher polydispersity of the suspension used here suppresses 
crystallization. On the other hand, the precise determination of the volume fractions at which crystallization occurs 
leads to a precise determination of the effective volume fractions (j)eff [Cheng et al. (2002)]. This fact has been used 
in our previous study using a narrowly distributed thermosensitive system [Crassous et al. (2006a, 2008a)]. In this 
way the zero-shear viscosity as well as the high-frequency and the high-shear viscosities could be obtained as the 
function of (j)eff- The effective volume fraction (j)f,ff may in principle be determined from the mass balance as lined 
out in previous studies [Crassous et al. (2006a, 2008a)]. However, the error in this determination turned out to be of 
appreciable magnitude. Hence, we followed the suggestion of Cheng et al. (2002) and compared the high-frequency 
viscosities of the narrowly distributed system studied previously [Crassous et al. (2008a)] with the respective values 
determined directly for the present system. Fig.|6]displays the resulting dependence of the high-frequency viscosities 
(hollow symbols) and the high-shear viscosities (filled symbols) on ^eff- The relative high frequency viscosity rj^ can 




be described by equation 4.22 in [Lionberger and Russel (1994)] with a reduced effective volume fraction of 1 1.95 % 
as opposed to the actual value of (jfeff, the volume transition neglegting a partial draining. The marked reduction of 
the volume fraction at high frequencies has already been detected in a previous study of a similar system [Deike et 
al. (2001)] and was explained by a partial draining of the thermosensitive network. The resulting correlation between 
the effective volume fractions and Tj^ of the naiTowly distributed systems can then be used to determine (j)eff of the 
present system. The values thus obtained are gathered in Table 1. Moreover, the inset of Fig.|6]demonstrates that there 
is a linear relationship between these ^^.ff and cRj^ as discussed above. A further check of the data is provided by a 

comparison of the high-shear viscosities rjL for both systems. Fig.|6]demonstrates that both sets of data agree within 
the given limits of error. This fact demonstrates that the present procedure leads to reliable values of (peff- 



Microscopic IMCT calculations 

While schematic models provide a handy description of the universal features close to the glass transition and the 
shear melting of glasses, the microscopic formulation of ITT-MCT actually aims to predict from first principles the 
nonlinear rheology starting with the underlying particle interactions; see the Appendix for a short summary. Because 
of the anisotropy of the structure under shear, it has not been possible yet to solve the ITT-MCT equations without 
approximations. Only solutions of schematic or isotropically averaged models are available. Yet, without shear the 
equations of quiescent MCT [Gotze 1991] can be solved readily and can be compared to the frequency dependent 
moduli within linear response regime. This has been done considering the monodisperse system in Ref. [Crassous 
et al. (2008a)], and after determination of the effective packing fraction we could perform identical calculations for 
the polydisperse system; see Ref. [Crassous et al. (2008a)] for details of the calculation which use the Percus-Yevick 
approximation to the structure factor of hard spheres. Here, of course the simplification to do monodisperse MCT 
calculations to compare with the measurements in the polydisperse system adds to the errors in the comparisons. 
The adjustable parameters, mainly the effective collective short time diffusion coefficient and the effective packing 
fraction setting the separation to the glass transition are reported in Table |2] Quite good agreement with the data 
can be achieved, with parameters that vary regularly with density and take quite reasonable values. Interestingly, the 
polydisperse system has larger moduli by a factor around two than the monodisperse system at the same effective 
packing fraction. 

TABLE 2. Experimental parameters (concentration c, temperature T and 
experimental effective volume fraction (pgjf with 0^ = 0.58. The parameters 
of the microscopic calculations (separation parameter e, volume fraction (j) 
with <pg = 0.5159, Cy and Ds/Dq, where Do = ksT / (GtztIsRh)), as shown 
in Fig. [5] are listed here. Furthermore, the power law exponents p from the 
dashed lines in the flowcurves of Fig.[5]are given in this table. 
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CONCLUSION 

The fluid-to-glass transition in colloidal suspensions leads to characteristic changes of the flow curves and concomi- 
tantly of G' and G" as the function of frequency. Here we demonstrated that both sets of data can be quantitatively 
modeled by a single generalized modulus g{t, 7). Moreover, g{t,Y) can be fully described by mode-coupling theory 
(MCT). Hydrodynamic interaction enters only through the high-shear and the high-frequency viscosities Tj,^ and Tj", 
respectively. MCT demonstrates that the latter quantities are not identical but closely related to each other (see the 
discussion of Eq.|7]). Small deviations between data and theory are only seen for G" in the immediate vicinity of the 



glass volume fraction and only at tiie lowest Peclet numbers Pea. This finding can be assigned to a very slow relaxation 
that is not captured by MCT. 
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Appendix 

The ITT approach, which generalizes the MCT of the glass transition to colloidal dispersions under stationary shear 
[Fuchs and Gates (2002); ibid. (2008)], consists of equations of motion for a density correlator <t>q(f ) which describes 
structural rearrangements, and approximated Green-Kubo formulae relating sttess relaxation to decay of density 
fluctuations. The approach was described in [Crassous et al. (2008a)], and the used equations can be summarized 
as follows. 

The density correlator obeys an equation of motion containing a retarded friction kernel which arises from the 
competition of particle caging and shear advection of fluctuations 

d,^^{t)+Y^{t) |<I>q(0 + jJ'^/f'mq(,-)(?-O5('^q(O|=0. (8) 

with the advected wavector q(?), whose component along the gradient direction varies with time, qy{t) = qy{0) — jtqx 
because of the affine component of the particle motion; while qx and q^, remain time-independent. The density 
correlator is normaUzed «I>q(f = 0) = 1 . Its initial decay is described by rq(/) = D^q^it) /•S^(f), with the collective 
short time diffusion coefficient, which may depend on hydrodynamic interactions, and thus differ from the Stokes- 
Einstein-Sutherland value. The equihbrium structure factor Sq, encodes the particle interactions and introduces the 
experimental control parameters like density and temperature. The generalized friction kernel Wq/(/), which is an 
autocorrelation function of fluctuating stresses, is approximated from the structural rearrangements captured in the 
density correlators 

'»q'(0 = ^1%^ Vp(OVq'kp(0) OkCO <I>p(0 , (9) 

where the vertices are given by the equihbrium structure factor 

Vqkp(0 = q(0 • (k(0«C<:(f) -|-p(0«Cp(,)) V+P ' (10) 

with n the particle density, and cj; the Omstein-Zemicke direct correlation function c;^ = (1 — 1/5;^) jn. Similarly, the 
generalized shear modulus is approximated assuming that stress relaxations can be computed from integrating the 
transient density correlations 

_ f klkyky{-t) S'kS'^^-,) 2 

^^''^>- 2 j {2nf kk{-t) Sl ^M-O(^)' 

where 5^ = dS^/dk. 

(y) 

The schematic -model described above follows from a simplification of these equations. It mimicks using 
Vi and V2 the increase of the vertices Vq'kp with sharpening structural correlations in Sq , which captures the cage 
effect causing a glass ttansition [Gotze (1991), Gotze and Sjogren (1992)]. It mimicks with 1/(1-1- {ft/Yc)^) the 
shear-driven decorrelation, which can be recognized from the vanishing of the memory kernel at long times under 
shear: mq/(/7^ 1) = and mq/(/) — m^(/)(^=*'' ^((yr)^) in the microscopic ITT equations [Fuchs and Cates 
(2002)]. Without shear, 7=0, Eqs. (8-11) become the famihar MCT equations which can readily be solved after 
discretization, as has been done in previous studies [Crassous et al. (2008a)]. The thin dashed lines in Fig. 5 labeled 
microscopic model result from these calculations, with parameters listed in Table 2. But under shear the equations 

(y) 

become anisotropic and require more severe approximations, Uke for example using the simple schematic -model. 

In detail, the equations were solved for a dispersion of monodisperse hard spheres with structure factor taken from 
the Percus-Yevick approximation. The wavevector integrals in the microscopic model Eqs. (8-12) were discretized 



according to Ref. [Crassous et al. (2008)] using M — 600 wavevectors chosen from k^i„ = 0.05 /Rh up to fcmax = 
59.95 /Rh with separation Ak = 0.1 /Rh- Time was discretized with initial step-width dt = 2 10^^/?^/Dj, which was 
doubled each time after 400 steps. Because the first principles calculations carry quantitative errors, which would mask 
the qualitative comparison, two adjustments of parameters were performed. The magnitude of the linear modulus is 
too small in theory by a factor Cy which is included in Table 2. The computed position of the glass transition is 
also somewhat to small, — 0.5159 in MCT relative to (pg = 0.58 in experiment. Therefore the relative separation, 
E = — ^g) / ^g was used as adjustable parameter quantifying the distance to the glass transition; the fitted values are 
also contained in Table 2. 
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